Random matrix averages and the impenetrable Bose gas in Dirichlet and Neumann 

boundary conditions 



P.J. Forrester, 1 ^ N.E. Frankel, 2 '0 and T.M. Garom 2 '0 

department of Mathematics and Statistics, University of Melbourne, Victoria 3010, Australia 
^School of Physics, University of Melbourne, Victoria 3010, Australia 
(Dated: February 4, 2008) 

The density matrix for the impenetrable Bose gas in Dirichlet and Neumann boundary conditions 
can be written in terms of (niLi I cos <^i — cos6l;|| cos 02 — cos#;|), where the average is with respect 
to the eigenvalue probability density function for random unitary matrices from the classical groups 
Sp(n) and + (2n) respectively. In the large n limit log-gas considerations imply that the average 
factorizes into the product of averages of the form (IIILi I cos0 — cos#i|). By changing variables 
this average in turn is a special case of the function of t obtained by averaging Y[?=i \t — xi\ 2q over 
the Jacobi unitary ensemble from random matrix theory. The latter task is accomplished by a 
^v^j , duality formula from the theory of Selberg correlation integrals, and the large n asymptotic form is 

obtained. The corresponding large n asymptotic form of the density matrix is used, via the exact 
solution of a particular integral equation, to compute the asymptotic form of the low lying effective 
single particle states and their occupations, which are proportional to yN. 
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PACS numbers: 03.75Hh, 05.30Jp, 02.30Rz 



> 1 I. INTRODUCTION 

| The probability density functions (p.d.f.'s ) 
t-H ■ 

- (2-) n 4sin2 ^) n 4(cos0 fe -cos^) 2 a) 

' ^ 71 ' ' 1=1 l<j<k<n 



Yl 4(cos6» fe -cos6i 3 ) 2 (2) 

l<j<k<n 



where < 0j < it (j = l,...n) occur in both random matrix theory and the quantum many body problem. In 
the former they are eigenvalue p.d.f.'s for classical groups with the Haar (uniform) measure - the group Sp(n) of 
n x n unitary matrices with real quaternion elements (which are themselves 2x2 matrices), and the group + (2n) 
k^j ' of 2n x 2n unitary matrices with real elements (real orthogonal matrices) and determinant equal to +1, for and 
@ respectively. A self contained derivation of these facts can be found in 0, Chapter 2]. In the latter they are 
the absolute value squared of the ground state wave function for n free fermions on the interval [0, 7r] with Dirichlet 
and Neumann boundary conditions respectively. As is similarly well known, and revised from first principles in our 
work 0, they are also the absolute value squared of the ground state wave function for n impenetrable bosons on 
the interval [0, tt] - in one-dimension the ground state wave function of the impenetrable Bose system is equal to the 
absolute value of the corresponding free Fermi system. 

In studies relating to both these seemingly disparate interpretations of the p.d.f.'s and J2J there is cause to 
investigate the function of (0, to, n) defined by averaging 



n 

1=1 



(cos0-cos6»;) m (3) 



with respect to these p.d.f.'s . In the random matrix interpretation this comes about in applications to L-function 
theory H, Q, Briefly, there are families of L-functions with special symmetries which are known to have their 
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non-trivial zeros well described by eigenvalues of random matrices from the classical group corresponding to that 
symmetry. For cos</> = 1 and small values of to the expected value of © can be computed with {6{\ corresponding to 
the zeros of particular families of L-functions, and it can also be computed - with m a general non-negative integer - 
for the random matrix ensembles. This then allows for both a test of the original hypothesis relating L- functions to 
random matrices, and provides specific conjectures for the statistical properties of the zeros of the L-function families. 
In the quantum many body interpretation the immediate interest is not in the average of but rather the average 

of 

n 

| COS <j)\ — COS 9i 1 1 COS (j)2 — COS 8l I . (4) 

1=1 

This gives the ground state density matrix of the corresponding impenetrable Bose gas system (if the absolute value 
signs are removed, the average gives the ground state density matrix for the free Fermi system) |2|. However the 
study of l@J leads back to the computation of J3J). Thus as noted in and for the problem of computing the 
asymptotic behavior of the ground state density matrix for the impenetrable bosons in the bulk and in an harmonic 
trap respectively, for large n and 4>\ and (\>i fixed the average (@J) is expected to factorize, and be proportional to the 
average of the product involving cf>i times the average of the product involving <p2- The latter are then the continuation 
in to of 10) from the even positive integers to the value to = 1 (this being a way to effectively study the average of 
the absolute value of ©)■ 

We remark that in the case of the density matrix for the impenetrable bosons in periodic boundary conditions the 
analogous task is to compute the average of 



m 



e »>i_ e «i|| e <Ai_ e «i| ( 5 ) 



with respect to the p.d.f. 



T| | e «*_e w '| a , (6) 

l<j<k<n 



where — 7r < 6j < 7T (j = 1, ...,n). In this case the asymptotic form for large n with 4>i and 4>2 fixed is a special case of 
known asymptotic forms for Toeplitz determinants with singular generating functions of the so called Fishcr-Hartwig 
type (for an extended discussion on this point and references to the relevant literature see Q)- The average of (0J 
with respect to Q or (J2J) can readily be written as a Hankel determinant, but there is no known analogue of the 
Fisher-Hartwig asymptotic form. 

In Scction[n]of this work we show how, for to an even positive integer, the average of © with respect to the p.d.f.'s 
(|T|) and (J2J), which is by definition an ri-dimcnsional integral, can be written as a m-dimensional integral. From the 
latter the large n asymptotic form of the average is deduced. In Section ITTT1 the result of Section [n] is used to deduce 
the large n asymptotic form of the average |0J with respect to the p.d.f.'s Q and J2J and thus of the ground state 
density matrix. We conclude in Section llVI by applying our result for the asymptotic form of the density matrix to 
the computation of the occupations of the effective single particle states. We also make some remarks in relation to 
the wider setting of our asymptotic analysis, in which we use Coulomb gas arguments to formulate an analogue of the 
Fishcr-Hartwig asymptotic form for a class of Jacobi unitary ensemble averages. 

II. DUALITY FORMULAS FOR MULTIPLE INTEGRALS AND ASYMPTOTIC ANALYSIS 

A. The duality formula 

The change of variables 

x . = i(cos^ + 1) (0 < Xj < 1, j = 1, ...n) (7) 
transforms and into the p.d.f.'s 

-A^) iimi^)) i/2 n (**-*j) a (8) 

'=1 l<i<fc<" 
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(=1 l<j<fc<n 



n! \2?r 

respectively. These p.d.f.'s in turn are special cases of the class of p.d.f.'s proportional to 



n^a-^^ n o»*-sj) 2 ( io ) 

1=1 l<j<k<n 

known in random matrix theory as the Jacobi unitary ensemble. Also, under the change of variables (JJJ the task of 
computing the average of © becomes the task of computing the average of 

n 

(=1 

with respect to (JHJ and © , or more generally with respect to (|l()|l . 

In fact there is an advantage in further generalizing the setting of the computation of the average of (|llfl and 
considering the class of multiple integrals known as Selberg correlation integrals, defined by 



* n / m 

- dx 1 ...dx n \\\x^(l-x l )^\\(t l 

6 J[o,i]» t\ \ 7=1 



S n ,m(^u ^2, A;t x , ...,t m ) dx 1 ...dx n | | ( x^il-xi) xi) 




(12) 



,7(2A)E„ 



Here 



C = S n (Xi,X 2 ,X)= d Xl ...dx n \[x^{l~ Xl ) x * J] l^"^| 2A ' (I 3 ) 

!=1 l<j<fc<ri 

known as the Selberg integral, is the normalization chosen so that the coefficient of ]X =1 1™ is unity, and the average 
over J(2A)E n refers to the p.d.f. 

1 n 

— — - — n^(i-xi) As n I**-*/* (1 4 ) 

6 n (Ax,A2,A)^ i</<fe<„ 

(the notation J(2A)E„ denotes the Jacobi -(2A) ensemble, which with A = 1 corresponds to the Jacobi unitary ensemble 
JH3J0- Setting t\ = ... = t m = t in lfT2^l gives the average of lfTT|l with respect to lfHjl. The advantage in studying 
(|12|l is that we can put to use the discovery of relating the Selberg correlation integrals to the theory of Jack 
polynomials (in the case A = 1 the Jack polynomials coincide with the Schur polynomials Q). In particular the 
Selberg correlation integrals were evaluated in terms of a generalization of the Gauss hvpergeometric function 2-F1 
based on the Jack polynomials. It was realized by one of the present authors 0, 0, El El El El that theory 
initiated in [S| could be further developed and used to express the average of © with respect to (114(1 and its limiting 
forms as the Laguerre - (2A) ensemble and the Gaussian - (2A) ensemble, as m - dimensional integrals. Because 
the role of n and m is effectively interchanged, these integration identities have been referred to as duality formulas 
One of their uses, as we will demonstrate in the case of the average of (fTTJl with respect to (fTT)|) . is in the 
computation of the large n asymptotics. 

The particular duality formula of interest to us is given explicitly in E3 • To state the result we must introduce the 
generalized circular ensemble, C/3Ejy, as the p.d.f. proportional to 

l[ \z k - Zj f, ( Zj = e^, -7T < 6 3 < 7r, j = 1, AO (15) 

l<j<k<N 
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With this notation, we read off from [15j eq. (3.41) that 



N \ I rn \ 

Y[zi m - m)/2 \i + Zl \^m +tz l rj oc (i - t)xi] N \ 

'=1 ' C^E N \i=l / J(4//3)E„ 



(16) 

Ai=2(Tja-Tn+l)/j8-l 
A 2 =2(7 7l +l)/ / a-l 

But we want to make (HlLi(^ ~ x l) m )j(2\)E ^he Q uan tity being transformed, so 1(16(1 requires manipulation. For this 
we write 

1 2 
ti-»l--, m^N, — = A, A = n (17) 



Noting that then 



i?i = i(A 2 + 1) - 1, »ft = i(Ai + l)+n-l, (18) 



multiplying both sides of 1)16(1 by i mJV and taking the complex conjugate of the left hand side of ((16(1 shows 

(f[(t - Xl) m ) = A m z ([(^-^)Al-")/ 2 |l + ^|[(Ax+A 2 +2)/A]+n-2 [i(1 + ^ _ (1Q) 

" =1 ' J(2A)E„ \i=l / C(2/A)E m 

where A, the proportionality constant, is independent of t. To specify A requires, in addition to the Selberg integral 
(pr5|> . the so called Morris integral 

— / d$!... d9 n \{z { r b)/2 \l + z l \ a+b [] \z k - Zj \ 2X . (20) 

^ ^ 1=1 l<j<k<n 

Then setting i = 1 in 1(19(1 shows that 

_ S n (Ai, A 2 + m, A) Af m (Q, 0, 1/A) 

5„(A 1 ,A 2 ,A) M m (7 72 , 7 7l ,l/A) 1 ' 

where rji and 772 are given by 1(18(1 . Both the Selberg integral and Morris integral have exact evaluations in terms of 
products of gamma functions (see e.g. Q). In the case A = 1 these read 

q( hU " TT 1 r(a + i + j)r(b + i + j)r(2 + J ) 

S n (a,b,l) = II — — — 7 

j=~o T(a + b + I + n + ]) 

G(n + 1 + a) G(n + 1 + b) G(n + 1 + a + b) + 



M„(o,6,l) = J] 



G(l + a) G(l + b) G(2n + 1 + a + b) 
1 T(a + b+l+j)T(2 + j) 



.L = l r(a+l + i)r(6 + l + i) 
_ G(n + l + a + fr) G(l + a) G(l + b) 

~ G(l + a + 6) G(n + l + a)G(n + l + 6) 1 + J ' 1 j 

where G(z) denotes the Barnes G-function, related to the gamma function by the functional equation 

G(z+l)=T(z)G(z). (24) 



B. Asymptotics 

Our interest is in the asymptotic form of the J(2A)E n average in ((19(1 in the case A = 1 and m even. The experience of 
our previous study , in which we studied the same product averaged over the Gaussian unitary ensemble (eigenvalue 
p.d.f. of complex Hermitian matrices with Gaussian entries) using a result known in the literature (see also [2(J), 
tells us the related quantity 

■^k,a 1 ,a 2 ((A, q)) 
^n+«,Ai,A 2 ((")0)) 
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where 

Z nMM {(X,q)) = X x ^(l-X) x ^ dXt... dX^X^il-X^X -Xtfi 

Jo Jo J=1 

x [] \Xk-Xj] 2 (26) 
l<j<k<n 

is better suited for the purpose. In addition to multiplying the JUE„ average in l|19[l by the i-dependent factor 
i Al " l / 2 (l — t) A2 ™ 1 / 2 , a key feature of (|25() is the normalization chosen so that in the interpretation of (|26() as the 
configuration integral for a log-potential Coulomb gas the normalization has the same total charge. 
It follows from (fH3)l . (|2~T|l . and noting 

-Zn + g,Ai,A 2 ((-, 0)) = 5 n + g (Ai, A2, 1), (27) 
that we have the duality formula for (|25|l with q = m/2, to even 



SV^Ai, A 2 4- to, 1) 



^uA 1 ,a 2 ((^,to/2)) = 

•^n+m/2,Ai,A 2 ((', 0)) S n +m/2 (Al , A 2 , l)Af m (Ai + 71, A2, 1) 



^Aim/2^ £^A 2 m/2 



y (Ai-A 2 -n)/2 
-/ 



A1 + A2 + U 



Ml + 2l) - 1]" 



1=1 



l<_7<A;<m 

SrJM, A2 + to, 1) 

(27ri) m 5„ +m/2 (Ai, A 2 , l)M m (Ai + n, A 2 , 1) 



4W, 



(28) 



where 



In® 



j\\raj1 



(l-t) 



A 2 m/2 



« m 

/ dx!...^™ TT(-^) Al+A2+ri (i - ^)" (A2+ " +m) (i - te;) r 



x (x k ~ Xj f 

l<j<fc<m 



(29) 



In l|29|) C is any simple closed contour starting at Xj — in the complex plane and encircling Xj = 1 anti-clockwise 
without crossing the interval Xj £ (0, 1). We obtain the second equality in (|28|l by writing the integrand in a form 
without absolute value signs, changing variables 



d6i = 



-dzj 



3 2mzj 3 



(30) 



then noticing the integrand is analytic except at Zj = 0, 1 and with a cut along the interval Zj € (0, 1). 

The large n, fixed t £ (0, 1) asymptotic analysis of an integral very similar to l|29l) has been detailed in [14|. and that 
analysis in turn follows the stationary phase analysis of a related multiple integral given in . Now, the n-dependent 
terms in the integrand of /„ (t) are 



x]{l - xj) "(1 - txj) n = exp[n(log£j - log(l - Xj) + log(l - tx ))]. 
As noted in |l4j , a simple calculation shows that the stationary point of the exponent occurs when 



l±i 



- t H-t) 



1/2 



X±. 



(31) 



(32) 



This suggests we deform the contours so that m/2 integration variables, [x\, x m / 2 ) say, pass through x+, and the 
remaining pass through x_. We must then expand the integrand in the neighborhood of these stationary points. 

Because we have made a definite choice of the to/2 variables, we must multiply by the combinatorial factor ^„^ 2 ^- 

From [1 1| we know that expanding the exponent in (|31|l to second order in (x± — Xj ) gives 



na 



exp[n(loga; J - log(l - x 3 ) + log(l - txj))] ~ exp[ — — (xj - x±) ] 



(33) 
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where 



{l-tx±Y 4 (l-x±) 2 ' 
Regarding the leading order expansion of the other terms in the integrand we have 

n (x k - X] f ~ (x+ - n k-^-) 2 n k-^o 2 

\<3<k<m l<j<k<m/2 m/2+l<j<k<m 



(34) 



(35) 



Hence 



II^)^^^ 1 _ ^)" (A2+to) ~ |a;+| (Al+A2)m |l - x+\- (X2+m ^ m . 

3 = 1 



J„(t)~ ( m ^ ) [i Al (l -<) A2 ] m/2 (x+ -x_) 2(m/2) |x + | (Al+A2)m |l-a; + r (A2+m)l 



n ^fc-^) 5 

l<j<k<m/2 



/oo />oo m/2 

dxi... / dx m/2 Y[ ex P 
-OO J — CO 

= [t Al (1 - t) Aa ] m/2 (a:+ - :r-) 2(m/2)2 M (Al+A2)m |l - x+\-^+ m ^ 

x |na|-( m / 2 ) 2 (Kn/ 2 ) 2 



where 



/oo ^oo 
dzi... / dx m/2 Y[ exp 
- oo J — oo 



II {xk-Xj) 2 

\<j<k<m/2 



From 



\ x +\ = \lj> li-^+l 



— — , |af+-af_| = 2W — — , |a| 



2t 3 / 2 



(l-t)V2 



so (|3*7ll simplifies to read 



/,-/) - (-l) m/2 ( J' ) 2^\-^ 2 [t{l-t)]- m2 l\V m/2 f 



(36) 



(37) 



(38) 



(39) 



(40) 



Furthermore, we recognize V m /2 as a limiting case of the Selberg integral known as the Mehta integral |2l| . which has 
the evaluation 



m/2-1 

V m/2 = (27r) m / 4 J] r(2 + j) = (2^r /4 G(m/2 + 2). 

3=0 



(41) 



Recalling l|28|l . our remaining task is to compute the asymptotic form of the combination of Selberg integrals and 
the Morris integral therein. According to (122[l and l|23ll, fo r this we require knowledge of the asymptotic expansion of 
the Barnes G-function. In fact Barnes himself showed [22J 



,'G(n + a + l)\ „ . a — b . i <r 

I,* l G \ n + b + 1 > ) n ~J> -*)n+ — log(2.) + ((a-b)n + 



logn + o(l). 



Using this we find 



SV^Ai, A 2 + 77i, 1) 



77. 



m 2 /2-m/2 



S , n+m/2(Ai,A 2 ,l)M m (Ai +77, A 2 ,l) G(m + 2) 



(42) 



(43) 
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Substituting (|41|l in 137(1 . then substituting the result together with (|43[l in (|28|) we obtain the sought asymptotic 
formula 

Z n+qMtX2 ((;0)) niG(2q+iy 1 U U V ' 

where we have set to/2 = q and use has been made of the functional equation l|24|) . We note that the right hand side 
of l|44|) is independent of the parameters Ai and A2 . 

A check on our workings to this stage is the special case q = 1- Then (|25|l coincides with the eigenvalue density in 
the JUE, p JUE (X), normalized so that its integral on [0, 1] is unity. Setting q = 1 in (|44|) we read off that 

p JUE (X)~-(X(l-*)r 1/2 , (45) 
which is indeed the known functional form (see e.g. Q]). 



III. ASYMPTOTIC FORM OF THE DENSITY MATRICES 



Consider the impenetrable Bose gas of TV + 1 particles confined to the interval [0,L] with Dirichlet boundary 
conditions. We know from p| that the ground state density matrix p^ +1 (x, y) is given by 



PN+i( x >y) = -^7 - ( sin ~]~) ( sIii t ! / j 



LJJo 



N 



dx N Yl 



1=1 



. 2 nx l 
sin — — - 

L 



TTX TTXl 
COS — COS 

L L 



cos 



ny nx l\ TT 1 nx k nx 3 

cos cos cos - 

L L I J-J- I L L 

l<j<k<N 



where 

C = 

Let us now change variables 



dx\. 



N+l 



dx N+ i J J sin 2 -j- J][ 

1=1 l<j<k<N+l 



TTXk TTXj 

cos — cos — j 1 - 



cos = 2Xj - 1 
L 



in both H46|) and (|47|l . and let us define 

f% +1 (X,Y) :=p° +1 (x,y) 



cos 7zx/L— 2X — 1 • 
cos Tty/L=2Y-1 



Then in terms of the generalization of (|26|l . 

Z nMM ((X, qi ),(Y, q2 )) = \ X -Y\ 2 ^X X ^(1-X) X ^Y X ^(1-Y) X ^ f dX x ...\ dX n 

Jo Jo 



n^ Ai (i-^) A2 i*-*;i 29i i^-^i 292 n \ x k-^ 

l<j<k<n 



1=1 



we have 



1 (*,n = TF ^[*(i-*)] 1/ W-n ^ 

\ X - y\ 1 Zjv+l,l/2,l/2((->°)> (.,0)) 



(46) 

(47) 

(48) 
(49) 



(50) 



(51) 



where p :— N / L. 

Similar considerations apply to the impenetrable Bose gas of N + 1 particles confined to the interval [0, L] with 
Neumann boundary conditions. Again from we know that the ground state density matrix p^ +1 (x,y) is given by 



Pn+i( x ,v) 



(N+l) 
C 



dx\. 







N 







1=1 



TTX TTXl 

cos — cos — — 



+ 

ny ttxi 

cos cos 

L L 



n'/IXfe VI Xj 
COS — COS - 

l<j<k<N 



(52) 
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where 



C = 



I dxi... I 
Jo Jo 



dx N+ i Yl 

l<j<k<N+l 



TTXk TTXj 

cos — cos — j 1 - 



Defining 



p% +1 (X,Y) :=(% +1 {x,y) 
and changing variables according to l|48[l in l|52[l and (|53|l shows 



COS 7TX / L — 2X — 1 
cos 7ry/L=2Y-l 



p% +1 (X,Y) = 



X-Y\V 



[X{1-X)]^[Y{1-Y)] 



1/4 ^,-i/2,-i/2((A-,l/2),(y,l/2)) 

^+l,-l/2,-l/2((-,0) ) (-,0)) 



(53) 



(54) 



(55) 



As already noticed in |(| Q , the log-gas interpretation of (|50ll allows us to predict that for large n it factorizes into a 
function of X and the same function of Y, which are themselves of the form (|26|) . Explicitly, we expect 



Zn > \i,\ 2 ((X,qi),{Y,q 2 )) Z nMM ((X,q 1 )) Z nM ,x 2 ((Y, q 2 )) 



Zn+ qi +q 2 ,\ u \ 2 ((',0), (•, 0)) Z n+qit Xi,X 2 ((•) 0)) Z n+q2) X u X 2 ((', 0)) 



(56) 



As in l|26[) . the key to choosing the correct normalizations is to balance the total charge in the log-gas interpretation. 
Setting qi = q 2 = q as required by 151|) and l|55|l it follows from (|56|) that 

Z nMM ((X q U(Yq)) ^ a^±l) {2n) - q+ A 2 [x{1 _ X)] -^ [Y{1 _ Y)] -^. (57 ) 

Zn+2 q ,X 1 ,X 2 ((-,0), (-,0)) \TTlG(2q+l) J 

Substituting this asymptotic form with q = 1/2 in (|51|l and (|55(l we obtain that for large N and fixed X, Y E (0, 1), 

G 4 (3/2) [X(l - AT)] 1/8 [r(l - y)]i/8 



I2N 



ix-y|V2 



(58) 



It is of interest to compare the asymptotic formula (|58J) against a numerical determination of say p® +1 (X, Y) or 
more conveniently p^ +1 (X, 1 — X). To compute the latter we write it as a random matrix average. Thus it follows 
from the various definitions that 



JV 



f% +1 (X, 1-X) = j^jX(1 - X) m (|4(1 - X) 4X,|) (\AX 4X t 



\/=i 



JUEj. 



(59) 



Ai = A 2 = l/2 



For each k = 1,2, ...M suppose we sample from JUEjv|ai=a 2 =i/2 obtaining the TV-tuple (x[ k \ X 2 , ...,X^'). Then 
the method of Monte Carlo integration tells us that 



M N 



p° +1 {X,l-X) = ^X(l-X)l^n(|4(l-X)-4Xf|)(|4X-4xW 



k=l 1=1 



O 



(60) 



Fortuitously, we have available a recently discovered [23j random three term recurrence which generates a polynomial, 
the zeros of which have the p.d.f. J(2A)E„. In the case of interest (Ai = X 2 = 1/2, A = 1) the recurrence states 



where with 



A (x) = 1 

A x {x) = x-BD[n+l/2,n+l/2] 

Aj(x) — {w 2 {x — 1) + w x)A 3 ^i{x) + wix(x — l)Aj^ 2 (x) (j = 2,...,n) 



aeGD[n + l- j + 1/2,1], 6eGD[j-l,l/2], c e GD[n + 1 - j + 1/2, 1] 



(61) 
(62) 
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L),MU(y -, Y\/ n O ( y 1 y\ 


0.025 


1.0958 


0.075 


1.0039 


0.125 


1.0363 


0.175 


1.0098 


0.225 


0.9439 






0.325 


0.9692 


0.375 


1.0338 


0.425 


0.9706 


0.475 


1.1309 



TABLE I: The ratio p°'^[ c (X, 1 -X)/p% +1 (X, 1 -X) where p^ c (X, 1-X) refers to the Monte Carlo expression while 
p^ + i(X, 1 - X) is the asymptotic form f5Sfy . We chose AT = 14, and evaluated JgSJ with Af = 5000. 



and rf := a + t + c we have 

w = 3 wi = 3, w 2 = 1 - wo - Wi. (63) 
a a 

Here BP [a, b] denotes the classical beta distribution, while GD[m, a] denotes the classical gamma distribution. The 
theory of 23] tells us that A n (x) has its zeros distributed according to JUE„| / \ 1 =a 2 =a- Implementing (|61|) for fixed n 
we thus computed the samples required in l|6(J|) for the Monte Carlo evaluation of pj$ +1 (X, 1 — X). Forming the ratio 
then with the asymptotic form l|58(l gave the data in Table [I] 



IV. PHYSICAL AND MATHEMATICAL IMPLICATIONS 

A. Ground state occupation of effective single particle states 

The ground state density matrix is the theoretical quantity which quantifies the condensation of a Bose system. 
Thus if we decompose the density matrix 

oo 

Pn(x,v) = ^2*j<i>j(x)<l>j(y), (64) 

3=0 

where the Xj , <f>j are the eigenvalues and normalized eigenfunctions in the eigenvalue problem 

l>yi r. g)(j) j (y)dy = Xj<f)j(x). ((>">) 

then by analogy with the free Fermi system in which (|64|) holds with Xj = 1 (j — 0, N — 1), Xj — (j > N), we 
see that the Xj have the physical interpretation as the occupation numbers of effective single particle states 4>j(x). 
For Bose-Einstein condensation to occur we must have Ao proportional to N . 

To study l|65|> in the case of the impenetrable Bose gas in Dirichlet or Neumann boundary conditions we restrict 
ourselves to large N where use can be made of the asymptotic form of the density matrices (|58H . Before mak- 
ing the substitution, recalling the definitions 149|) and l|54|) we must first change variables in (|65|l and redefine the 
eigenfunctions so that 

COSny/L = 2Y - 1, <t>j{Y) = <t>j{y)\ co a-ny/L=2Y-U </>j( X ) = <l>j(x)\ cos7TX / L= 2X-l- (66) 

Doing this we obtain, for large N, the integral equation 

~2 n J \x~Yjn ^ (y Vy(i-y) = M y (67) 

It follows immediately that 



X j oc VN. 



(68) 
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As noted in 0] this conclusion requires that j be fixed - for j N we expect \j oc (N/j) 4 in keeping with the 
corresponding result in periodic boundary conditions, since in this regime the boundary conditions are not expected 
to play a role. 
Setting 

G 4 (3/2 ) 

and rearranging, (|67|l reads 



Xj = p > VN\ 3 (69) 



" =A JT ^2) . (70) 



o IX-FI 1 ^ [y(i-y)]3/s 3 [X(1-X)]V8- 

Remarkably the effective single particle ground state </>o(X), and the corresponding scaled occupation number Ao can 
be computed exactly from (|70|l . To see this requires knowledge of a piece of integral equation theory presented in 
Porter and Stirling |2_^ . The relevant theory tells us that the solution of the integral equation 

0(<) dt = l, v<l (71) 



o V-A v 



is 



7T 

Setting v — 1/2, it follows immediately that 



ifcos^Wl-zr- 1 )/ 2 . (72) 



b {X) = ^=[X{l-X)] l l\ Ao = Try/2 (73) 



and so 



satisfies (|71|l . where the normalization ^4 is determined by the requirement that 

A = —5(3/4,3/4) (75) 

where B(a,b) denotes the beta function. Substituting the exact evaluation of A in (|69[) shows that in the large N 
limit 

A = G 4 (3/2)\/iV = 1.3069V77. (76) 

To compute the higher order single particle states and their corresponding occupations we make the ansatz 

<f>j{X) tx <fo(X) Pj (X) (77) 

where Pj(X) is a polynomial of degree j. Now {(j)j(X)} can always be chosen to be orthogonal (note that the measure 
is dX/\JX{l — X) on [0, 1]) so recalling lf?3")l we require 

rM^) dX = 3+k . (78) 
Jo (X(1-X))V4 ^ > 

Up to normalization, the unique polynomials with this property are the particular Gegenbauer polynomials 

Pj (X) = C^(2X-l), (79) 
which we note are proportional to the particular Jacobi polynomials P- 1 ^ 4: 1 ^ 4 (2X — 1). Normalizing l|77|l with the 



substitution (|79")) as in (|7U) shows 
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Now substituting (|80() in (|70|l and setting X = 1 we obtain a definite integral for Xj which can be found in [25| , giving 
us the evaluation 

Xj = ^ F(j+ , 1/2) (81) 

and hence 

Xj = G 4 (3/2) r(j + 1 { 2) Vn. (82) 

To arrive at l|8UI) we have made the ansatz (|77H . In fact a different approach can be taken to the problem, in which 
it is shown that an integral operator following from l|70l) commutes with the differential operator determining the 

polynomials {C 1 ^ i {2X — l)} :)= o ! i.2,...-This is done in Appendix 1X1 

Finally, we note that substituting (|58|) . I|8UI) and l|81|l in (|64H gives the following interesting identity 

J / x-^, . .,.,-,1/4, ov 1^^1/4/ 



y_ y |i/ 2 = V^ 2(1/4) EC? + V±)C l /\2X - 1)C) / \2Y 1). (83) 

I j=0 



B. Generalized Fisher-Hartwig type asymptotics 

One viewpoint of our asymptotic analysis of multiple integrals of the form (|50|) is that we are studying asymptotic 
problems of the Fisher-Hartwig class. Let us recall that the latter refers literally to Toeplitz determinants with both 
zeros and jump discontinuities is its generating function, 

oo 

D N [e a W] := detta^-]^!,...^, e a(e) = E a p e^ e (84) 



P— — OQ 



where 



with 



a{9) = g(6) - i E M 71 " ~ (9 - r )]mod27r + E a r log |2 - 2 cos((9 - r )| (85) 



r— 1 r— 1 



7(<?) = £ 2p e ^> E bll<7p| 2 < oo. (86) 



p= — oo p— — oo 



Thus g(0) is a regular term, while at <j) r (r = 1, R) there is a jump discontinuity of strength b r and a zero of order 
a r . To see the relationship with l|5U|l we set b r = (r = 1, R) thus eliminating the jump discontinuities, and recall 
the general formula relating a Toeplitz determinant to a multiple integral, 



JSk _ e iSj |2 



' •'° ^° 2 = 1 l<j<k<N 

1 /.27T /.27T N / R x 

= -M ^n^' ni eie, - e ^i 2a n 



2 = 1 V=l / l^X/c^N 



(87) 



Fisher and Hartwig [2(| conjectured that in the case (|87|l . 

Z>jv[e a(e) ] ~ e ^ e £-i a > sJV i; (88) 
where is independent of N . This was subsequently proved, and it was furthermore shown 

E = e T.T=l k 9k9-k TJ e -M9(<M-3o) TJ | e *^* _ g#i |-2a fe a.j -Q G 2 (l + Or) 
r=l l<j<k<R r=l ~^ ^ ar ^ 
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(see for example the monograph [27j and references therein). 

As noted by one of the present authors some years ago [2^|, it is straight forward to reproduce the structure of 
lf£%|l using the analogous log-gas argument to that used here in the analysis of lf5T)|l . Now the (normalized) multiple 
integral corresponding to l|87|l which relates to 1501 is 

HnMM[e hix) Utl\yr-^] (Qn) 

f7 m ' ^ ' 



-1 .1 n 

, ' ^° 1=1 l<j<k<n 



\ Xk Xj | 



(91) 



where h(x) is analytic on (0, 1). As our final issue, we would like to extend the log-gas argument used in the analysis 
of l|5U|) to predict the large n asymptotic form of (|90|) . 

From the log-gas perspective, the natural quantity to analyze is 



II \ yk - yi 

l<j<k<r 

where for m a non-negative integer 



a q , qk HnMMle h(x) T\Li\yr-x\ 2qr } 



„1 „\ m 

H m [h}:= dxi... dx m J[e h ^x^{l-xi) x * ]J \x k - Xj f 

l<j<k<m 



1=1 



Analogous to (|56|l we expect for large n (|92|l to factorize as 



HnMM^ nti \Vr ~ s| 2 *1 A - qM H nM [\y r -x\ 2 ^} 



H 



n +E ^ iqM2 W l{x) ] 



n 



H n+qr: \ u \ 2 [l] 



(92) 



(93) 



(94) 



where the second expression is motivated by inspection of the known results i|88|) and (|89|l for (|87|l . 
Thus we expect 



H n ,x u \ 2 [^] 



-2 qjqk _ 



H, 



l<j<k<r 

qr h{y r ) Hn,\iM [\Vr ~ x \ 2Qr 
H n+qrt x u \ 2 [l] 



+E?=i 9r,Ai,A 2 



>] 



>< 



But according to (fP^) 



^A^Rlfr-af*] 1 G 2 (q r + 1) 



(2n)-^(y r (l-y r ))-^ 2 . 



(95) 



(96) 



ffn+fr.A^p.] irlr G(2q r + 1) 

Also, for the first ratio in (|94|l we have available both rigorous results [2^, [3(J as well as log-gas type arguments 
which together tell us that 

H n+QMM [e h ^] 



-ffn,Ai,A 2 [l] 



exp 



x exp 



x exp 



n + Q + (Ai+A 2 )/2 f 1 h(x) 



[x(l - x)] 1 ' 2 



dx 



Ai + A 2 



(MO) + &(!)) 



dx- 



h{x) 



_4tt= J [ X (l- X )]i/2J 
Substituting ll9l)|) and (|9"?|l in (jMJ gives the analogue of (|55j) , 



(97) 



-^n^i.AsIl] 



exp 



x exp 



E?=i<?r + (Ai+A 2 )/2 Z" 1 



J2(~q r + q 2 r )log2n K 



o [a:(l-x)]Va 



(98) 
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where 



K = II I?/* — 2/j -|- 2 ^^ e -C A i+ A 2)[/ l (0) + ft(l)]/4 e -i:« =1 
l<j<k<R 



x exp 

R 



1 

4^ 



dx- 



h(x) 
x(l -x)] 1 ' 2 Jo 



dy- 



x-y 



x 



1 G 2 (g,- + 1) 
LA nir G(2q r + 1) 



n 



(99) 



C. Concluding remarks 

In our first paper on the impenetrable Bose gas Q we set ourselves the goal of providing the leading asymptotic 
form of the density matrix for the impenetrable Bose gas in a harmonic trap and in Dirichlet and Neumann boundary 
conditions. It was noted in [2| that for the impenetrable Bose gas in periodic boundary conditions, the Fisher-Hartwig 
formula gave the asymptotic form 

r , , G 4 (3/2) / 7T \ 1/2 , s 

In , it was shown that for the harmonic well 

(2N)^ P % +1 (V2NX, V2NY) ~ N ^Hm (1 - X *?'*(l-pV* (m) 

while in the present paper, after changing variables according to (|49(l . the asymptotic form of the density matrix 
is shown to have the leading asymptotic form (|58J) . As a consequence of the scaling properties of these asymptotic 
forms, the occupation number of the low lying effective single particle states are all proportional to yN, but with a 
proportionality constant dependent on the particular system. 

To obtain the asymptotic forms we have used a combination of exact analysis, made possible by the theory of 
Selberg correlation integrals, and physical reasoning based on log-gas analogies. Taking this argument to its logical 
conclusion leads to a conjectured exact asymptotic formula, given by l|98|l and (|99|l for a Jacobi weight analogue of 
the Fisher-Hartwig formula. 



APPENDIX A: PROOF THAT {<t> {] (X)C) /A {2X - l)} 3 =o,i,2,... ARE SOLUTIONS OF THE INTEGRAL 

EQUATION (1701) 

The assertion that the <j>j{X) given by l|77|l are solutions of 17U|) is equivalent to stating that the Gegenbauer 
polynomials are eigenfunctions of the integral operator 

*vw := L w^W^Wr^ (A1) 

where for convenience we are working on the interval — 1 < £ < 1. In this Appendix we prove that K commutes with 
the differential operator, L, which determines the Gegenbauer polynomials 

^-P -<">£- 14 (A2) 



with 



LCy i (0 = -3ti+W- /4 (0. (A3) 
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We begin by identifying if[C^ 4 (£)] with the following finite sum of hypergeometric functions 



K[cy\o] = a, E 



4(-i) fe (j + l/2) fe /1 + A 1/4 _ (I . 35 1 + C 



fc=0 



fc!(3/4)» 



2*1 



^ e 



H)*(7 + V2)fc A-£\ 1/4 i-i/l ,3 5 1-$ 



fc=0 



fc!(3/4)» 



where 



r(3/4) r(j + 1/2) 

r(5/4) j! 



(A4) 



(A5) 



To derive (|A4(I we break up the interval of integration in K[Cj (£)] into two regions, thus removing the need for the 
modulus, which gives 



1/4 



v-1/4 „l/4. 



i-iy 4 /4 /'#r 1/4 (i-^r 1/2 (i-^)- 1/4 ^ /4 (i-2^) 

Jo 

where we have defined 



2± 



1±C 



(A6) 



(A7) 



and in obtaining the second equality we have changed integration variables from ip to 

(1 — ?/>)/(l — an d (1 — + m the first and second integrals respectively. Substituting the following known 

25] expansion for the Gegenbauer polynomials 



(-l) J T(j + 1/2) ' (-j) fc (j + l/2) fe (1 + $ 



3W n 



E 

k=0 



fc!(3/4) fe 



into (|A6|) . and using the standard integral representation of the 2-F1 function 

# LwILm = 5(3/4, 1/2)2^ - - A, -; -;z 



(1 - z^)- fe+1 / 4 



4' 4' 



(A8) 



(A9) 



then results in (|A4|) . 

We can now utilise known hypergeometric identities to facilitate the operation of the differential operator L on the 



expression for -^[C ? 1//4 (£)] given by (|A4|) . We note that the structure of (| A4|) is of the form 



K[cy 4 m = + i-iy^sjiz-), 



where 



a (z) ._ y (-i)fc(j + 1/2)* .1/4 pi /I , 

,() -~ fc tj *!(3/4) fc * 



3 5 



'4'4 ;Z 1 ■ 



In terms of both the variables z = z + and z = z_ the differential operator L has the form 

d 2 3. 



d 
dz 



Utilising the identity (3S 

d n 
dz n 



z° F 



«i, a p 

Pi, -,Pq 



(S-n+l) n z s ~ n p+1 F q+ 



(5 + 1, ai, a p 
5 + 1 -n, pi, ...,p g 



(A10) 



(All) 



(A12) 



(A13) 
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we find that 



4 ' 4' 4' J z V 16 / V 4 4 4 



(A15) 



^(i-4B = 4 1/4fc ^(z-4i ; * 

-r 1/4fc2Fi (i-^li^ 

_ z i/4 fc (I _ ^ 2Fl Q _ (fc _ 1); 3. 5. ^ (Alg) 
where the equalities in (|A15fl and (|A16|I respectively follow from the two particular contiguity relations [32] 

■h)^-^ Fi (\- k 'h-b z ) = hw-")*-*un(z- k >vv* 

- %*Fi(i-k,r,r,*), (ait) 



2 * 1 V 4 4 4 
and 

-J^l (J - *, | 5 ^) - g-*,^.,) *) ^ (i - (* - 1), »; ,) . (A18) 

Hence from (|A16|I we obtain 



fc=0 



fc!(3/4) fc V 4 / V 4 4 ' 4 



i -2 



Changing summation index in the second sum in l|A19|l from k to k — 1 and simplifying we deduce finally that 

rcM , , /0 s ( _ j)*C7 + V 2 )fc 1/4 77 Z 1 ; 3 5 

LS,(z) = - j(3 + l/2)g fc , (3/4)fc ,/ 2 F^--fc,--- 

= + l/2)^(z), (A20) 

which then implies that 

^C 1/4 (£) = -i(j + l/2)^Cj /4 (0 (A21) 
= KLC)!\i) (j = 0,1,2,...) (A22) 

and so 

[if, L] = 0. (A23) 
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